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The purpose of this study, which is based on photographic recordings showing the 
phase shift of two signals — i.e., GBR transmitted from Rugby (England) on 10 kilocycles 
per second and NBA transmitted from Balboa (Panama) on 18 kilocycles per second, both 
received at the Neuchatel Cantonal Observatory (Neuchatel, Switzerland), is twofold: 

(1) To investigate the mode of propagation of the GBR and NBA signals. 

(2) To study the relation between the time of sunrise, respectively sunset a1 various 
ionospheric reflect ion points and the times at which phase fluctuations appear on the re- 
cordings. The author then generalizes the notion of the times of sunrise, respectively 
sunset by introducing the closely feinted concept of the zenithal distance Z of the sun at 
the reflection points considered. Following results published in literature, reflection point 
altitudes were assumed to be about 70 kilometers during the day. Results were such that: 

(a) For the GBR signal; only a one-hop mode is available, night reflection altitudes 
varying between 88 and 91 kilometers on individual recordings, angles of incidence <£ on 
the ground between 7°36' and 10°25'. 

(b) For the XB A signal; a five-hop mode is available, night reflection altitudes varying 
between 80 and 84 kilometers on individual recordings, angles of incidence </> on the ground 
between 0°27 / and 1°14'. 

(c) At sunrise, respectively, sunset, computed zenithal distances for one and the same 
reflection point at times identical with singularities appearing on successive recordings 
show a st riking analogy. 

(d) During one and the same sunrise or sunset, the zenithal distances computed suc- 
cessively for various reflect ion points and related to singularities read on a same recording 
present also a striking analogy. 

(e) The time of onset of ionizing radiation at all nighl reflection points seems to be 
of major importance for both the propagations of the GBR and NBA signals. During 
sunset, the altitude of the day reflection point which was stabilized at around 70 kilometers 
increases as soon as the zenithal distance of the sun exceeds 90°. At sunrise, on the other 
hand, the altitude of the reflection point stabilizes at around 70 kilometers, when the zenith- 
al distances of the sun reach or go below 90°. 

(f) In the case of the NBA signal a phase fluctuation already occurs at a zenithal 
distance of about 103°, especially at sunrise. At that moment, the distance between the 
reflection point and the layer formed by the ionizing radiations of the sun is about 100 
kilometers. 

(g) At sunrise, the curves showing the energy of the received signals display Hie follow- 
ing features : 

For a one-hop mode (Rugby), a momentary strong absorption when the reflection point 
altitude reaches 82 kilometers; in the case of several ionospheric reflections (Balboa), a 
succession of absorption lines corresponding to the successive diminishing of the altitude 
of the ionospheric reflection points. 

(h) At sunset, the interpretation of the energy is more delicate: 

For a one-hop mode (Rugby), a momentary increase occurs in the energy of the signal 
before the night level of reflection is reached; in the case of a five-hop mode 4 (Balboa), the 
interpretation of the absorption curve is difficult because five ionospheric reflection points 
change their altitude and the resulting phase fluctuations become entangled. 



1. Introduction 

A simple geometrical model explaining the transi- 
tion from night-to-day propagation resulted from 
a previous experimental study [Rieker, I960]. The 
model assumed longwave propagation in successive 
hops in the waveguide formed by the earth and the 
lower edge of the ionospheric layer. Implicit in this 
was the assumption that the angles of incidence on 
the ground are constant for a definite wavepath. 



Atmospherics emanating from the same azimuth 
had been picked up on the loop aerial of a narrow- 
sector direction finder pointed toward the source; 
these receptions were analyzed from recordings. 

This study confirmed a view previously held, that 
reflection-point altitudes are higher at night than 
during daylight hours. Moreover, it became possible 
to distinguish mixed propagation paths from all- 
night propagation or all-day propagation. It was 
also found that angles of incidence on the ground are 
very small, varying between and 2 to 3°. 
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In the meantime, theoretical studies, especially 
Wait's fundamental "Diffraction Theory for LF 
Skywave Propagation" [Wait, 1961; Wait and 
Conda, 1961], indicate that the angles of incidence 
on the ground should be in excess of about 5° for a 
27 kc/s wave frequency if geometrical optics are 
strictly valid. Despite violations of this condition, 
geometrical optical arguments are often used. From 
a qualitative viewpoint, the final conclusion should 
be valid. For a more refined treatment, rigorous 
mode theory should be used [Wait, 1962]. For 
example, it was found in the study mentioned above 
[Rieker, 1960] that 27 kc/s waves emanating from 



distant thunderstorms (sferics) are received at angles 
of incidence on the ground less than 2 or 3°, the mode 
of propagation being the simplest possible, as pre- 
viously suggested by Wait [1959]. 

Obviously, the method used in investigating dis- 
tant thunderstorm emissions cannot be applied 
at sunset. Narrow-sector direction finder recordings 
[Rieker, 1960, (figs. 1 and 2)] show that these sources, 
when recorded at night, usually cannot be located 
during daytime, at least as far as more distant thun- 
derstorms are concerned. This results in major 
difficulties when an attempt is made during the day 
to determine the azimuth of areas of thundery 
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Figure 1. A diagram of the photographic method used at the 
observatory of Neuchdtel for the registration of the phase fluc- 
tuations of the GBR signal (Rugby) and NBA signal (Balboa) . 

Accuracy of the atomic clock: 1 part in 10". Frequency difference between 
oscillator in Balboa and atomic clock in Neuchatel always smaller than 2 parts 
in 10 10 (i.e., 17 ^sec per 24 hr). 
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Figure 2. Phase fluctuations of the GBR signal at sunset on the Sept. 10, 1960. 
Discontinuities A, B, C, F, and G are recorded in table 12. 
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activity to be recorded the following night. Further- 
more, it is impossible to check selected azimuths 
without operating a narrow-sector direction finder 
under normal conditions. 

Here we propose to reinvestigate the problem of 
longwave propagation by another method: the 
analysis of phase fluctuations of a signal (Rugby 
GBR or alternatively Balboa NBA) with respect to 
a local frequency standard (atomic clock at Neu- 
chatel). 

Scrutiny of recordings kindly supplied by Dr. J. 
Bonanomi, Director of Neuchatel Observatory, re- 
vealed the possibility of some conclusions about 
longwave propagation at sunrise and sunset. 

2. Instrumentation and Recording 

As a check against the frequency of the Neu- 
chatel standard (atomic) clock, the Neuchatel Ob- 
servatory (47°00 / — 06°57'E) receives signals from 
Station GBR, Rugby, England (52°22' N— 01°11' W) 
and Station NBA, Balboa, Panama (09°03'N— 
79°39 / W). These frequencies are intercompared 
(fig. 1) by means of a filmed record of the oscilloscope 
trace, which is triggered by the Neuchatel reference 
signal. 

The cathode tube screen shows a fixed sinusoid 
when the phase difference is constant. At sunrise 
and sunset, the length of the signal path is modified 
by the change in altitude of the ionospheric reflec- 
tions points. A phase lag at sunset is indicated by 
a shifting of the sinusoid, with an equal and opposite 
shift at sunrise. Filming of this sinusoid shift at 
two or three points (1, 2, 3) of identical ordinate 
c results in nearly parallel traces, as shown in 
figure 1. 

Let T=- be the period of the sinusoid, where 
v 

v is the frequency of the received wave, o)=2ttv 

is the gyratory frequency, and c is the distance from 

the slit to the abscissa. The following equation 

expresses the sinusoid on the oscilloscope screen: 
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T f is expressed in microseconds and is related to the 
period by the equation 

The square of amplitude 
A 2 =- 



cos^ 



T 



(3) 



reveals the behavior of the received wave energy. 

When S designates the total path traveled by the 

wave after completion of n successive hops, the 



length difference AS between the day and night 
paths is derived easily from the equation giving the 
phase: 

/ T' y 



which can be written 



5 = 27rf, 

A 



and finally 



AS = vA(t l + 



('•■ 



v 



(4) 



v being the velocity of propagation of the signal, 
i.e., 3.10 5 km/sec= velocity of light. 

The accuracy of the atomic clock is 1 part in 
10 11 . The frequency difference between the oscil- 
lator in Balboa and the atomic clock in Neuchatel 
never exceeds 2 parts in 10 10 (i.e., 17 jusec per 24 hr). 

Four recordings of the phase fluctuations of the 
GBR signal (figs. 2, 3, 4, and 5), and five recordings 
of the phase shifts of the NBA signal (figs. 6, 7, 
8, 9, and 10) were studied. These recordings are 
among the first made at the Neuchatel Observatory 
during initial trials, which explains the two different 
scales used for time registration. 
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11.9.60.-5unrise : Neuchatel - Rugby 

Figure 3. Phase fluctuations of the GBR signal at sunrise on 

the Sept. 11, 1960. 

Discontinuities I, K, L, and M arc recorded in table 10. 



3. Computation of the Length Difference 
AS (km) Between the Day and Night 
Paths 

3.1. Rugby 

For Rugby, the period of the sinusoidal wave is 

r=r / +T"=62,5/isec. 

The following values were obtained when comput- 
ing the differences AS=vtj (4) between the fixed 
points A— G, figures 2 and 4, for sunset, and I—M, 
figures 3 and 5, for sunrise: 



Figure 2, sunset 
Figure 4, sunset 
Figure 3, sunrise 
Figure 5, sunrise 



A£ 7 =31 jusec, 

A^/=28 jiisec, 
A£ 7 =28 jusec, 
A£ 7 =26 jusec, 



AS=9,3 km 

AS=8,4 km 
AS=8,4 km 
AS=7,8 km. 
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2A161.-Sunset: 
Neuch^tel- Rugby 

Figure 4. Phase fluctuations of the GBR signal at sunset on 
the Jan. ££, 1961. 

Discontinuities A, B, F, and G are recorded in table 12. The square of the 
amplitude of the signal represented below has been calculated from the time 
difference T", resp. T"hi/xsec, (see also fig. 1). 
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25.1.61.-Sunrise: 
Neuchatel-Rugby 

Phase fluctuations of the GBR signal at sunrise on 
the Jan. 25, 1961. 

Discontinues I, K, and M are recorded in table 10. The square of the ampli- 
tude of the signal represented below has been calculated from the time difference 
T'in Msec, (see also fig. 1). 



Figure 5. 



Computation^of^the path length differences as- 
sumes a signal propagation velocity equal to the 
velocity of light. 

3.2. Balboa 

For Balboa, the period of the sinusoidal wave is 

T=T'+T" =55,55 /xsec. 

Since these recordings (figs. 7 and 9) were inter- 
rupted before sunrise had occurred over the complete 
path Neuchatel-Balboa, they cannot be considered 
for the time being. The following values were 
obtained when computing the differences AS=vt T (4) 
between the fixed points A x and i^* 5 , figure 6, B x 
and F 5 , figure 8, and ^4*! and G 5 , respectively H x and 
P 5 , figure 10: 

Figure 6, sunset A£ 7 = 59 Aisec, AS=17,7 km 

Figures, sunset A/ 7 =49 ptsec, AS=14,7km 

Figure 10, sunset -f sunrise A£ 7 = 71 jusec, A>S=21,3km. 

4. Determination of the Mode of Propaga- 
tion 

4.1. Altitudes of the Day Reflection Layer 

In order to determine the mode of propagation of 
the signal, the altitudes of either the day or night 
reflection points must be known. Since neither can 
be found in this study, reference was made to results 
of previous research found in literature indicating an 
altitude of day reflection layer of about 70 km 
[Budden, Ratcliffe, and Wilkes, 1939; Bracewell and 
Kruse, 1951; Nertney, 1953; Gregory, 1956; Eevellio, 
1956; Wait, 1959]. 

This being assumed, nocturnal reflection altitude 
computations were based on the values of AS 
deduced from the readings of A£ 7 on the recordings (4). 

4.2. Mode of the Rugby GBR Signal Propagation 

Table 1 provides AS values as computed for a 
one-hop mode, figure 11, and for a two-hop mode. 
For the particular case of only one ionospheric 
reflection two neighboring: reflection altitudes were 
selected, i?=88 km and 77=90 km. 

Table 1. Neuchdtel — Rugby 
(Geocentrical angle: 0/=7°3O') 
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Figure 6. Phase fluctuations of the NBA signal at sunset on the Sept. 7, 1960. 
Discontinuities arc recorded in table 16. 
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8.9.60.-Sunrise : Neuchatel- Balboa 

Figure 7. Phase fluctuations of the NBA signal at sunrise 
on the Sept. 8, 1960. 

I >iscontinuities are recorded in table 14. The sunrise-effect continues after the 
Interruption of the record at 9 UT. 
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Figure 8. Phase fluctuations of the NBA signal at sunset on the Sept. 11, 1960. 
Discontinuities are recorded in table 16. 
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12.9.60.-Sunrise : Neuchatel - Balboa 

Figure 9. — Phase fluctuations of the NBA signal at sunrise on the Sept. 12, I960. 
Discontinuities are recorded in table 14. The sunrise-effect continues after the interruption of the record at 9 UT. 
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Figure 10. — Phase fluctuations of the NBA signal at sunset on the April 22, 1961 and at 

sunrise on the April 23, 196L 

Discontinuities are recorded in table 16 and table 14. The inverse square of the amplitude of the signal represented 
below has been calculated from time difference T", resp. T" in /xsec, (see also fig. 1). 



It can be ascertained that if a one-hop mode 
occurs, the difference AjS=7,81 km, respectively 
8,76 km is comparable with the differences deduced 
from the recordings (see sec. 3.1) and are between 
7,8 km and 9,3 km. Angles of incidence <t> on the 
ground meet Wait's requirement [Wait, 1961]. 

A two-hop mode would result in a difference in 
the path traveled by the signal about 3,5 times 
larger. This possibility can therefore be excluded 
as well as that one of two-hops during the day and 
one-hop during the night or inversely. 

Ultimately therefore, the mode of propagation 
of the GBK signal received at Neuchatel is also 
only a one-hop mode, night reflection altitude 
varying between 88 and 91 km for a day reflection 
altitude assumed to be at 70 km. 

4.3. Mode of the Balboa NBA Signal Propagation 

a. Method 

The characteristic elements relating to the path 
Neucbatel-Balboa are shown in table 2. These ele- 
ments were computed for 5, 6, and 7 hops and also 
to a certain extent for 8 hops. The altitude of the 
day reflection point was taken as 70 km while two 



night reflection altitudes— i.e., 90 and 84 km— were 
considered. The column headed "5 ionospheric re- 
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Figure 11. Model of the propagation between Rugby and 
Neuchatel in the case of a one-hop skywave (night propagation 
1, day propagation 2). 
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Table 2. Neuchdtel — Balboa 
(Geocentrical angle: 0/=81°O5') 
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on the 
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9190,1 





5°53' 


084=5 O O4 / O4 // 


0*J= 5 o 04'04 /, 


573, 21 


9171, 4 







16070 =81°05' 


4°35' 


7 o=5°O4'O4" 




570, 72 


9131, 5 


58,6 


4°35 / 


7 o=5°O4'O4" 




570, 72 


9131, 5 


39,9 
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flections" for instance shows first the relation between 
geocentrical angles for the paths 1 and 6 (night and 
day propagation) as well as 2, 3, 4, and 5 relating to 
so-called intermediate propagations, figure 13. Then 
angles of incidence </> on the ground, geocentrical 
angles 6 70 , d 91 respectively 84 , as well as the geo- 
centrical angle for the last night hop with respect to 
Neuchatel, E | 0&=1» 2, ... 5), are presented, s is 
the length of the path from one reflection point on 
the ground to the next following on the lower edge 
of the ionosphere. S=2s the total length of each 
peculiar path 1, 2, 3, 4, 5, and 6, figure 13, AS being 
the difference in the length of each one of them with 
the length of path 1. 

Figure 12 illustrates the graphical method used 
to determine the above elements. The functions 
8== f(fi) an d s = f(4>) are given for three different 
altitudes. A third diagram (0 7O — 084, 90) shows the 
constructed curves <£ 7 o, 90 and </> ro , 8 4, locus of points 
such that for the corresponding <£ 70 and </> 8 4.9 angles 
of incidence on the ground $ 70 and <£ 84) 90 are equal 
(see graphical construction of the angles of incidence 
<£ 7 o,9o, respectively 7O 84 =0° and cfoo.oo, respectively 
^70,84=4°). The straight lines numbered 1 to 15 
are the relation between geocentrical angles with 
regard to the geocentrical angle Neuchatel-Balboa. 
Each bundle of straight lines represents one mode 
of propagation. 

The following example is an illustration of the 
method, figure 12: the straight line 9 intersects the 
curve 7O , 90 at K in the diagram (0 7O — 084, 90). K is 
first related to the abscissa 84 . 90 where 9O =7 o 38 / 45" 
is found; afterward .990=860,0 km is read on the 
curve £T=90 km in the diagram (s—6). This value 
on the corresponding curve H=90 km in the diagram 
(5 — 0) yields the angle of incidence on the ground 
0r=2°O2 / on the abscissa. The path length s 70 = 
738,0 km is found by intersecting the projection line 
with the curve H=70 km in the diagram (*,0). Re- 
turning to diagram (s,0), 7O is found to be 6°34'45". 

This method has the advantage of being quick. 
It can furthermore be applied to all propagation 
modes, provides of course, that the straight line 
equations are adapted first to the geocentrical angles 
as given by the locations of both transmitter and 
receiver. The graphical method just described is, of 
course, not as accurate as the mathematical method. 
Night and day path lengths for each mode 5, 6, 7, 
and 8 were accurately calculated. Global differences 
between night and day paths are therefore known 
with sufficient accuracy, while partial differences be- 
tween intermediate path lengths are far less accurate 
since they were found with the above graphical 
method. 

b. Determination of the Mode 

It was established in section 3.2 that length dif- 
ferences of the path traveled by the NBA signal 
varied between 14,7 km and 21,3 km for day and 
night propagation. 

If it is now assumed that night reflection occurs 
at an altitude of 90 km and day reflection at 70 km 
according to table 2, a difference of 31,5 km for 5 



hops, 39,1 km for 6 hops, 47,2 km for 7 hops, and 
58,6 km for 8 hops is found. For a reflection altitude 
of only 84 km during the night, these differences are 
respectively 21,6 km, 26,8 km, 32,8 km, and 39,9 
km. It may be noted that on the registration of 
April 22-23, 1961, figure 10, the difference AS=21,3 
km is similar to that for 5 hops with a reflection alti- 
tude of 84 km. Inferred differences for the examples 
of September 7, 1960, figure 6; and September 11, 
1960, figure 8; yield reflection altitudes of about 82 
and 80 km for the same mode. 

Consequently, it seems fairly well established that 
mode 5 is the actual propagation mode of the NBA 
signal between Balboa and Neuchatel. 

The problem has to be considered more closely, 
however. Assume, for instance, a propagation of 
mode n during the night (n hops) and a propagation 
of mode 71+ 1 during the day (n-\- 1 hops) or inversely. 
It becomes immediately obvious that the eventuality 
of having a mode 71+ 1 during darkness and n during 
daylight must be eliminated, the resulting differences 
of path lengths AS being far higher than those de- 
duced from the registrations. The other eventuality, 
n ionospheric hops during darkness and n-\-I hops 
during daylight, is more logical and looks tempting 
at first sight. 

A scrutiny of the global differences of the traveled 
paths for H=90 km and i?=84 km (as extracted 
from table 2) shows the following: 

5 night reflections} A&> = 91 16,5-9099,3 = 17,2 km; 

6 day reflections/ AS 84 =9106,6-9099,3 = 7,3 km 

a . , . n ,. 1 A# 90 = 9138,4-9114,8 = 23,6 km; 

6 night reflections fa 4=9 126,0-9H4,8=il,2 

7 day reflections 1 8 ~ 

AS 90 =9162,0-9131,5 = 30,5 km; 
AS 84 =9147,6-9131,5 = 16,1 
km. 



7 night reflections 

8 dav reflections 



The differences thus calculated correspond closely 
to those deduced from the recordings, especially 
for 17,2 and 23,6 km. 

However, transition from one mode to another 
can only take place in a precise order. 

If an intermediate propagation of mode n (n 
ionospheric reflections) is considered, n D being the 
number of day reflection points and n N the number 
of night reflection points, the following relation 
obviously applies: 

n=n D +n N . 
The mode of next higher order would be 

7l'=7l+l. 

Transition from mode n to n' must follow the law 

n'=(n D +l)+n N 

where n D +l=?7^, number of day reflections for 
mode n f . 
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Figure 12. Graphical method to find the angle of incidence <f>, the distance s between the iono- 
spheric reflection point and the next ground reflection point, and the corresponding geocentrical 
angle 6 for a determined mode of propagation that is represented by a bundle of straight lines. 



Legend: equations of the straight lines for 



5 ionospheric reflections 

1:2070+80*4 90= 81°05' 
2:407O+S084'9O=81°O5' 
3:607O+4084*9O=81°O5' 
4:80 7 o+2084 , 9O=81°O5 / 



G ionospheric reflections 
5:2070+10084 90=81°05' 
6:4070+8084 90 =81°05' 
7:60 7 o+6084'9O =81°05' 
8:8070+4084*90 = 81°05' 

9:10070+2084*90 =81°05' 



7 ionospheric reflections 
10:2070+12084 90=81°05' 
11:407O+1O084'9O=81°O5' 
12:6070+8084 90 = 81°05' 
13: 8070+6984*90 =81°05' 
14: 10070+4084*90 = 81°05' 
15:12070+2084*90 = 81°05' 
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IONOSPHERE 



E 1 E 2 E 3 




Balboa 



Neuchatel 



, IONOSPHERE 

D 6 ,^ 1 
X \ -J 




-6 f = 81°05' 

Figure 13. — Geometrical model of the longwave propagation about sunrise in the case of a 5-hop 
sky wave for the distance Balboa-Neuchdtel (night propagation 1, intermediate propagation 
2. 3. 4- 5, day propagation 6). 



Figure 13 shows this behavior more explicitly. 
Path 3 for instance (intermediate propagation 3) 
involves 2 day {D\ and Dl) and 3 night reflection 
points (Ef, El, and El), so that 







n D =2 








n N =3 








n = 5. 




The transition to n' = 6 iono 


sph 


such that 










™D 


= (n D +l) 


= 3 




nir 


=n N 


-3 




n' 


=n+l 


-6 



because the decrease of the altitude of the last night 
reflection point Ef occurs simultaneously with the 
addition of a new reflection point to the remaining. 

In table 2, path 3 discussed above is shown by 
(2) (No. 2 in fig. 12) while (7) (No. 7 in fig. 12) 
expresses the transition. 

What now are the effects of a transition, in accord- 
ance with the law indicated above from mode n to 
n' on the behavior of S, the length of the path trav- 
eled by the signal? Table 2 shows that whatever 
the altitude of the night reflection chosen (H=90 or 
84 km) and regardless of the timing of the transi- 
tion from one mode n to the next higher mode 
n'=n-\-l, the length of the path S increases. In- 
versely, if a night as well as a day propagation are 
of the same mode (e.g., paths 1, 2, 3, 4, 5, and 6, 
fig. 13, and table 2) (5 hops), the length of the path 
is always shorter than that immediately preceding. 

If, then, we assume a change of mode occurring 
together with sunrise over the path Neuchatel- 
Balboa, a sudden increase in path length should 
take place at the moment of change of the mode. 

The recordings, figures 6 to 10, however, show 
that the shift Ah of the sinusoid always occurs in 
the same direction at sunrise, respectively at sunset; 



therefore the respective variations AS of the path 
length traveled by the signal are similar. While at 
sunrise the length S of the path traveled decreases 
more or less evenly showing no tendency to in- 
crease at any time, at sunset, on the other hand, the 
length increases with no sign of decrease at any time. 

It can be concluded accordingly, that the mode of 
NBA signal propagation is the same both at day and 
at night. Mode five is therefore found to be appli- 
cable, altitudes of reflection points being about 70 
km during daylight and between 80 and 84 km at 
night according to the recordings available, figures 
6/8, and 10. ^ 

Angles of incidence on the ground for a five-hop 
mode in table 2, for a day altitude of 70 and 84 km 
at night, then, are 

O°27 / <0<1°14 / . 



4.4. Consideration on the Night Reflection 
Altitudes of the GBR and NBA Signals 

The night reflection altitudes of the GBR signal 
were found to be higher than those of the NBA 
signal. 

Since both frequencies (16 kc/s for Rugby, 18 
kc/s for Balboa) are quite close to each other, night 
reflection altitudes would be equal for identical 
geometrical propagation conditions. This being not 
the case, the inequality 

B E \\ (Balboa) >0* | (Rugby) 

must be borne in mind. 

This implies that angles of incidence on the lower 
edge of the ionospheric layer are also unequal: the 
GBR signal will penetrate further into the ionosphere 
and reflect at a higher altitude than the NBA signal. 
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5. Introduction to the Detailed Study of the 
Geometrical Model of Longwave Propa- 
gation at Sunrise and Sunset 

The recordings for Rugby reveal a sudden trans- 
lation of the sinusoid on the oscilloscope screen 
from / to M at sunrise, figures 3 and 5. On the 
other hand, sunrise recordings for Balboa are differ- 
ent in this respect. On April 23, 1961 for instance, 
5 successive shifts numbered 1 to 5 in figure 10 are 
shown. 

The question as to the relation between the sudden 
shifts of the sinusoid and the ionospheric reflection 
points now arises, since the altitude decrease of each 
of them must result in a shift of the sinusoid on the 
oscilloscope screen. A correlation between the times 
of sunrise at the reflection points and the times at 
which the shifts occur could be searched. However, 
computing a sunrise time at a reflection point is 
equivalent to a computation of the zenithal distance 
Z of the sun at this point. This new notion will 
furthermore be used to fix similar singularities of 
possibly identical origin on various recordings as will 
be shown later. It can be assumed that similar 
singularities related to the same reflection point, but 
taken successively on different recordings must by 
necessity have equal zenithal distances. When com- 
puted under identical conditions for several reflection 
points in a multihop mode propagation, similar 
zenithal distances will again be found both for sun- 
rise and sunset. The zenithal distance would there- 
fore appear to be independent of the declination of 
the sun and the position of the reflection point. 
1 This is not unlikely, since the rate of production of 
electrons in the ionosphere is known to vary with 
the cosine of the zenithal distance of the sun. 

This is the approach to the detailed study of the 
proposed model, figure 13, which will also clarify 
the model related to sunset, figure 14. 



6. Geometrical Considerations 
6.1. Rugby 

It is established in section 4.2 that the mode of the 
GBR signal propagation is only a one-hop mode. 
Let iV be the midpoint of the arc of great circle 
joining Neuchatel and Rugby along the ground, 
figure 11, D the day and E the night reflection 
points. Starting from the coordinates of Rugby 
(52°22 / N— 01°11' W) and Neuchatel (47°00' N— 
06°57' E), it is found that 

N (49°16' N— 03°16' E), 

D (49°16' N— 03°16' E— 70 km), 

E (49°16' N— 03°16' E— 88 km), 

with d f =7°30', the geocentrical angle between the 
receiver and the transmitter 



6 e — &D — @N 



11=1° AK' 



3°45', 



the azimuth of Rugby from Neuchatel being 

A Z =31S°37'. 

The recordings, figures 2, 3, 4, and 5, had shown an 
altitude of the night reflection layer varying between 
88 and 91 km (see sec. 4.2). This altitude was 
selected at 88 km for subsequent computations. 

6.2. Balboa 

a. Sunrise 

The intermediate propagation (paths 2, 3, 4, and 5, 
fig. 13) requires the reflection points on the iono- 
spheric layer contained in table 3 (azimuth 

^ z =273°44 r ). 



IONOSPHERE 




Figuke 14. Geometrical model of the longwave propagation 
about sunset in the case of a 5-hop skywave for the distance 
Balboa- Neuchatel (day propagation 1, intermediate propaga- 
tion 2. 3. 4- 5, night propagation 6) . 
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The geocentrical angles B B \ . . . Qe\ formed by 
each of the points E\ . . . El tabulated between the 
two diagonals of table 3 with the receiver (Neu- 
chatel) were calculated in table 2 on the basis of the 
geocentrical angles 6 70 and B M found with the 
graphical method already described, figure 12. 



Table 3 



Path 



Reflection point 




In a previous study [Rieker 1960, sec. 2.2 and 2.13], 
a simple method of computation of these geocen- 
trical angles was evolved, their general expression 
beino; 



_(2m-l) a[n-(m-l)} 

\ m-l/£ V D m/L2m—1 



(2m-l)+a, 



(6 n n+l\ 



2n 



(5) 



where a is defined by the relation 

(B^rn} — (Or*™) = <X, 
\ E \ JL \ D m/L ' 

for our case a^0 o 45' , n being the number of iono- 
spheric reflections (n=5 in our case, see sec. 4.3b), 
d f =81°05' the geocentrical angle between the trans- 
mitter (Balboa) and the receiver (Neuchatel); m an 
integer defined by the inequality l<^m<n, the 
index L relating the expressions to sunrise. 

When the azimuth and geocentrical angle of each 
characteristical point tabulated between the diago- 
nals of table 3 are known, the coordinates A and <p 
of each point can be computed easily; these are 
shown in tables 4 and 5. 

The recordings (figs. 6, 8, 10) yielded altitudes 
of night reflection varying between 80 and 84 km 
(see sec. 4.3b); a value of 84 km was selected for 
further computations. 



Table 4 



Reflection 

point 
77=84 km 


6 


<P 


X 


E) 


8°07' 


46°55'N 


4°56 / W 


E\ 


23 16 


43 40 


26 03 


E\ 


39 03 


36 42 


44 38 


E\ 


55 23 


27 00 


60 11 


E\ 


72 22 


15 26 


73 34 



Table 5 



Reflection 

point 
ft =70 km 


e 

y4*=273°44' 


<P 


X 


D\ 


7°31 / 


46°59'N 


4°05'W 


D\ 


22 59 


43 45 


25 41 


I>1 


39 03 


36 42 


44 38 


D\ 


55 43 


26 47 


60 28 


D\ 


72 58 


14 52 


73 53 



b. Sunset 

It was shown in section 4.3 that a five-hop mode 
must be applied to the Balboa-Neuchatel transmis- 
sions. This involves consideration of an intermedi- 
ate propagation 2, 3, 4, and 5 for the transition from 
day to night propagation similar to the case of sun- 
rise considered previously, assuming angles of in- 
cidence on the gTound for each particular path 1, 2, 
3, 4, 5, or 6, figure 14 to be identical. 

Table 6 shows all sunset reflection points. By 
reasons of symmetry the geocentrical angles between 
Neuchatel and the reflection points tabulated be- 
tween the two diagonals of table 6 can be derived 
from (5): 

2 m— 3 

where ; as previously in the case of sunrise, 

n being the number of ionospheric reflections and m 
an integer, the index C relating the expressions to 
sunset. 

The azimuth ^4 Z =273 44' and the geocentrical 
angle of the characteristical reflection points with 
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respect to Neuchatel being known, the coordinates 
A and (p of each of them can be computed easily. 
These are shown in tables 7 and 8. 

Table 6 



Reflection point 




Table 7 



Reflection 

point 
h =70 km 



A,=273°44' 


<p 


X 


D\ 


8°07' 


46°55'N 


4°56'W 


D\ 


25 22 


42 55 


28 45 


D\ 


42 03 


35 04 


47 43 


D\ 


58 06 


25 14 


62 29 


D\ 


73 34 


14 35 


74 27 



Table 8 



Reflection 

point 
11= 84 km 


e 

^ r =273°44' 


<P 


X 


E\ 


8°43' 


46°53'N 


5°50'W 


El 


25 42 


42 40 


29 06 


E\ 


42 03 


35 04 


47 43 


E\ 


57 49 


25 25 


62 15 


El 


72 58 


14 52 


73 53 



7. Further Developments 

7.1. Methods 

a. Introduction 

Prior to the detailed study of individual recordings, 
a summary of the three variables used in their 
interpretation and the methods of their computation 
follows: 

The zenithal distance Z of the sun, 

the times of rising and setting of the ionizing 

radiation of the sun, 
the distance d between one ionospheric reflection 

point and the layer formed by the ionizing 

radiations of the sun. 



b. Determination of the Zenithal Distance of the Sun 

The following computation method was used to 
determine the zenithal distance Z of the sun for a 
given time T and an ionospheric reflection point E 
respectively D of coordinates (A, <p, h), figure 15; 
from 

tj R 
cos U= lf¥Ii 



where 



and 



U= angle of depression 
i?=radius of the earth =6371 km 



— sin £7= sin 8 sin <p+cos 8 cos <p cos (£ + 0> (7) 

the value of 

(to+t) LtC 

is calculated for sunrise L or sunset C at point E or D 
of coordinates (X, <p, h). 

8 is the declination of the sun, t the hour angle of 
the sun at sunrise or sunset, t the time period between 
either sunrise or sunset and the time considered. 

For a known time difference 

L-T or C~T 

the following is derived easily: 
at sunrise 



(t +f) T =(to+t) L +(L-T) 



at sunset 



(to+t) T =(t +t) c -(C-T) 
The zenithal distance is given by 

cos Z=sin 8 sin <p +cos 8 cos <p cos (2o+£) a 



(8) 



c. Determination of the Time of Rising or Setting of the Ion- 
izing Radiation cf the Sun 

Sunrise or sunset at any given point (X, <p, h) 
near the surface of the earth is associated normally 
to the appearance or disappearance of the solar 
disk in the visible spectrum. However, the main 
ionizing agent of the high atmosphere is the ultra- 
violet radiation of the sun of below about 3000 A; 
it is logical, therefore, to assume a photoionic effect 
of this radiation on the air-components of the high 
atmosphere, when investigating the relationship 
existing between the times of sunrise or sunset at 
points along the path of a signal and the times at 
which phase fluctuations of this signal appear on 
the recordings. It is also known that this ultra- 
violet radiation is completely absorbed by the ozone 
layer. In other words, this means that this ultra- 
violet radiation will reach the point (X, <p, h) after 
the sun has risen and will leave it before the sun has 
set in the visible spectrum. The problem is, there- 
fore, reduced to a computation of the times of 
sunrise and sunset in the visible spectrum referring 
to a sphere with a radius of R+A r reaching the 



66&696— 63- 
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6 ^ 5 \2» 



§ angular 
A6'\*" sunradius 




1: transmitter 
2:receiver 



SUN 



Figure 15. One-hep skywave: zenithal distance Z and dis 
tance d /rom ^e reflection point E or D to the surface formed 
by the ionizing rays which are tangential to the ozone layer. 



sphere of influence of the ozone layer. The intro- 
duction of the notion of the rising and setting of the 
ionizing radiation of the sun is thus justified. 

A simple computation method for the time differ- 
ence results from a known altitude A r of the sphere 
of influence of the ozone layer above the surface of 
the earth. Starting from the equations 

where Ui and U 2 are the angles of depression, R the 
radius of the earth =6371 km and h the altitude of 
the point under study. 

Equation (7) then yields the difference of the 
times 

|(*o+*i)-(*o+* 2 )|=A. 

The time at which the ionizing radiation rises at 
an ionospheric reflection point thus becomes 

i+A, 

the time of their setting becoming 
C-A. 



The altitude A r being difficult to determine on 
the recordings under study, reference was made to a 
previous study [Rieker, 1960, sec. 1.3] where it was 
found to be 28 km; this value was adopted for 
further use in the present study. 

d. Computation of the Distance d From One Ionospheric Reflec- 
tion Point to the Layer Formed by the Ionizing Radiations of 
the Sun 

Some phase fluctuations of the GBR or NBA 
signal occur before the ionizing radiations of the sun 
have risen at the corresponding reflection point 
(X, cp, h). In this case the distance d from the re- 
flection point to the layer formed by the ionizing 
radiations has been calculated, figure 15. 

When, at sunrise, the phase fluctuation of the 
signal occurs before the ionizing radiations of the 
sun reach the reflection point considered, the rate 
of electron density must have been modified before 
the said radiations reached the particular reflection 
point. 

The layer formed by the ionizing radiations 
issuing from the sun taken as a punctiform light 
source assumes, in the first approximation, the shape 
of a cylinder tangential to the ozonosphere and 
therefore to the sphere R-\-A r , figure 15. 

It is easily recognized that for a ionospheric 
reflection point E respectively D, d is defined by the 
equation 

d= (R+A r ) - (R+h) sin Z. (9) 

Z is the zenithal distance of the sun at the day or 
night reflection point considered; it should not be 
mistaken for the zenithal distance of the sun at 
Neuchatel at the same moment, which would have 
an altogether different value. 

If, for curiosity's sake, the solar disk in the visible 
spectrum is assumed to have a diameter of 16", the 
distance d would have to be corrected by a value I 
determined by the expression (see fig. 15): 



l=(R+h) cos (ir-Z)tg 16'. 



(10) 



Assuming that: 



d 



100 km 
A r =28 km: altitude of the sphere of influence of the 

ozone layer [Rieker, 1960, sec. 1.3] 
5=6371 km: radius of the earth 
h=S8 km: altitude of the reflection point (see 
sec. 6.1), 

it is found that 

Z=102°47' 
Z=6,6 km. 

Conversely, the variation of Z resulting from 
d=6,6 km is Z^15'. 

As will be shown in the next paragraph, this error 
can be neglected, the inaccuracies resulting from the 
readins-s taken on the recordings being of the same 
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order of magnitude, sometimes even more important. 
The sun can therefore be taken as a point light source. 
It might be added that a strict reasoning would 
involve the dimensions of the sun in the ultraviolet. 

7.2. Interpretation of the Recordings 

a. Introduction 

The Rugby recordings will be studied first, 
figures 2 to 5, followed by the more complex recording 
from Balboa, figures 6 to 10. Recordings at sunrise 
being easier to interpret, they will be examined first. 

Mention has already been made of the fact that 
these recordings were not made with the specific 
view of studying long wave propagation. This 
explains the interruption in the recordings of Sep- 
tember 8 and 12, 1960, figures 7 and 9, before the 
sun had risen over the whole path Neuchatel-Balboa. 
Furthermore, since these recordings were part of the 
initial trials conducted at the Neuch&tel Observatory, 
two different speeds of advance of the 35 mm film 
were used; on the original high speed recordings, 
figures 2, 3, 6, 7, 8, and 9, 1 mm represents 2 min, 
while on the low speed recordings, figures 4, 5, and 
10, 1 mm corresponds to 8, 5 min. A reading 
error of 1 min results in an error for Z of AZ^30' 
in the first case, while on the low speed recordings, 
a reading error of 4 min results in an error of 
AZ^1° to 2°. 

The "Tables Crepusulaires" [Lugeon 1934, 1957] 
were used to compute the times of sunrise and sunset 
at a point (X, <p, h). 

It should be mentioned finally that the January 
24 and 25, 1961 recordings, figures 4 and 5, were used 
earlier in a preliminary study [Rieker, 1961]. 

b. Interpretation of the Recordings of the Rugby GBR Signal 

Recordings of September 11, 1960, Figure 3 and 
January 25, 1961, Figure 5 at Sunrise. In spite of 
two different time scales having been used (abscissa) , 
both recordings show a sudden initial phase shift 
at / followed by a clearly marked stabilization at K 
in figure 3 and M in figure 5. From point M onward, 
no significant shifting of the phase occurs on both 
recordings. Table 9 shows the time of rising of the 
sun and of the ionizing radiation at the characteristic 
points N, D, and E defined in section 6.1, figure 11, 
as well as the corresponding zenithal distances Z 
defined in section 7.1b, figure 15. 

Table 9 



Place 


Time of sunrise 


Time of incidence of the first 
ionizing sunray 




11. 9. 60 
UT 


25. 1. 61 
UT 


Z 


11. 9. 60 
UT 


25.1.61 
UT 


Z 


Neuchatcl 


5h09«» 

4 28 
4 21 


7hHm 

7 35 
6 37 
6 31 


90°00' 
90 00 

98 27 

99 28 








N 








D 


4h4()m 

4 32 


6*>50™ 

6 42 


«.)(•) XV 


E 


97 49 







Table 10 comprises the results as computed for 
the limes /, K, L, M indicated on figures 3 and 5. 



The zenithal distances Z of the sun at the times 
/, K, L, M were calculated at point N, figure 11 for 
both recordings. At a given time, K for instance, 
the zenithal distance is identical for all points on 
a vertical line from point N. The distance d 
between the night reflection point E (49°16 / N — 
03°16' E— 88 km) and the layer formed by the 
ionizing radiations was calculated for the time 
/ (see sec. 7.1d, (9)). 

Table 10 

Sunrise 



Abscissa 
in the 


Reflection 
point 


11. 9. 60 


25. 1. 60 


figures 3 
and 5 


UT 


Z 


d 


UT 


Z 


d 


J. 


E 


4h29m 

4 40 

5 03 
5 31 


98°19' 
96 33 
92 59 
88 27 


km 
8 



6^39™ 
6 57 


98°19' 
95 30 


km 

8 


K. 


D 

D 




L.. _.. 




M 


i).__ _ 




7 37 


89 42 









The following considerations are based on table 10: 

(1) On both recordings (September 11, 1960 and 
January 25, 1961), the initial shift of the phase 
occurs at a time corresponding to the time of rising 
of the ionizing radiation of the sun at point E 
(49°16' N— 03°16' E— 88 km). 

(2) Both recordings show a singularity at the 
times K coinciding with a zenithal distance of the 
sun of about 96°. 

Plotting A 2 ((3) sec. 2) below the January 25, 1961, 
recording as a time function calculated by repeated 
measurements of 7 7/ Gusee) an absorption curve 
showing the variation of the energy of the GBR 
signal is obtained. A minimum appears in this 
curve at time K which occurs at a reflection point 
altitude of about 82 km. Further singularities are 
noticeable in the absorption curve at times / and K. 

(3) The singularity indicated at time L, figure 3, 
is not visible on the recording of figure 5. This fact 
is probably due to the slow speed of advance of the 
film. 

(4) From time M onward i.e., when the ionizing 
radiation reaches the day reflection point D prac- 
tically horizontally, the altitude of the point of 
the ionospheric reflection has nearly reached the 
lowest level. The received energy, on the other 
hand, has slightly increased again (see absorption 
curve, fig. 5). 

Recordings of September 10, 1960, figure 2 and 
January 2\, 1961, figure 4 at sunset Phase fluc- 
tuations of the GBR signal are less clearly dis- 
cernible at sunset than at sunrise. The recording 
of September 10, 1960, for instance, shows a very 
gradual phase shift beginning at time A, whereas 
the beginning of the phase shift is slightly clearer 
on the recording for January 24, 1961. The phase 
shift, furthermore, extends over a much longer 
period of time at sunset than at sunrise and generally 
appears to be more complex. Beginning at time G, 
however, the phase shift becomes very small. The 
increase of the altitude of the night reflection point 
is probably terminated at this time G. 
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As previously for sunrise, times of setting of the 
sun and of the ionizing radiation as well as the 
zenithal distances Z of the sun as computed for the 
characteristical points N, D, E defined in section 
6.1 are presented in table 11. 







T 


A.BLE 11 








Place 


Time of sunset 


Time of incidence of the last 
ionizing sunray 




10.9.60 
UT 


24.1.61 
UT 


Z 


10.9.60 
UT 


24.1.61 
UT 


Z 


Neucbatel 


17h49m 

18 06 

18 59 

19 06 


16hl7m 

16 22 

17 19 
17 26 


90°00' 
90 00 

98 27 

99 28 








N 








D 
E 


18h4 7 m 

18 56 


17W"' 
17 15 


96°33' 
97 49 



Table 12 comprises the zenithal distances of the 
sun and the distances d from point E (49°16 / N— 
03°16'E — 88 km) to the layer formed by the ionizing 
radiations at times A, B, C, F, G indicating singu- 
larities apparent on the recordings. 

The following considerations are based on table 
12: 

(1) At times A, figures 2 and 4, the zenithal 
distance is about 90° for both recordings. From 
that time onward the phase shift becomes apparent, 
i.e., the reflection point altitude begins to increase. 

(2) When the zenithal distance approaches 94°, 
the phase shift of the GBR signal becomes slower. 

(3) A singularity seems to be significant at time 
C on the September 10, 1960 recording, figure 2, 
this particular time C being the time of setting of 
the ionizing radiation of the sun at point E 
(49°16'N— 03°16'E— 88 km) . This singularity does 
not appear on the recording of January 24, 1961, 
figure 4, or is at best only barely recognizable. 

(4) When the zenithal distance reaches about 102°, 
a new singularity appears in the phase shift (time 
F) . The distance from point E to the layer formed 
by the ionizing radiations then is about 100 km. 

(5) Beginning at times G, the altitude increase of 
the night reflection point seems to be completed. 
It should be pointed out here that the zenithal 
distance of the sun and the distance from point E 
to the layer formed by the ionizing radiations on 
September 10, 1960, and January 24, 1961, are 
too dissimilar to postulate the existence of a relation 
between the two. 

(6) By plotting A 2 ((3) sec. 2) below the recording 
of figure 4 as a time function, deduced from repeated 
measurements of T', respectively T" (/zsec), it would 
appear that critical points on the absorption curve 
are associated with times A, B, F. At time B a 
momentary increase of reflection capacity is notice- 
able. The reflection point would then have reached 
an approximate altitude of 73 km. Near time F 
absorption would have diminished considerably. 

Discussion. A first conclusion attained by com- 
parison of the results obtained at sunrise and sunset 
is that the altitude of the ionospheric reflection 
point is around 70 km as long as the zenithal distance 
of the sun lies between and 90°. Zenithal distances 
at time A at sunset and time M at sunrise are about 
90° according to tables 10 and 12. 



Table 12 



Sunset 



Abscissa 
in the 


Reflection 
point 


10.9.60 


24.1.61 


figures 
2 and 4 


UT 


Z 


d 


UT 


Z 


d 


A 
B 
C 


D 
D 
E 
E 
E 


18^06™ 
18 27 

18 58 

19 30 
22 48 


90°00' 
93 26 
98 19 
103 06 
125 09 


km 

______ 

108 
1118 


16^16^ 
16 55 


89°29' 
94 45 


km 


F 

O 


17 41 

18 50 


101 41 
112 42 


74 
440 



At sunrise, the phase shift seems to be triggered 
by the rising of ionizing radiations at point E 
(49°16'N— 03°16'N— 88 km), the beginning of the 
phase shift occurring at times /. 

At sunset, the layer formed by the ionizing radia- 
tions again seems to be of same importance, the 
distance of the night reflection point E from this 
layer being about 100 km at times F while the 
ionizing radiations disappear at time C, figure 2. 

The significance of the singularities noticed at the 
intermediate times K (sunrise) and B, F, G (sunset) 
becomes evident when related to the behavior of 
the energy A 2 , figures 4 and 5. 

When several ionospheric reflections (multihop 
mode) of the signal occur, the phase shifts which 
are completed in a relatively short time at sunrise, 
will be easily dissociated and analyzed. At sunset, 
on the other hand, the phase shift continues over a 
longer time, so that the interpretation of the re- 
cordings is more difficult; each peculiar phase shift 
produced by a successive increase of altitude of the 
reflection points being entangled with the next one. 

c. Interpretations of the Recordings of the Balboa NBA Signal 

Recordings of September 8, 1960, Figure 7; September 12, 1960, Figure 
9; and April 23, 1961, Figure 10 at Sunrise 

Unfortunately, the September 8, 1960, figure 7; 
and September 12, 1960, figure 9 recordings were 
terminated before sunrise had occurred over the 
whole path Neuchatel-Balboa. The recording of 
April 23, 1961, figure 10, on the other hand, is com- 
pleted. Two different speeds of advance of the 
film were also used. 

A first examination of the recordings of the NBA 
signal at sunrise shows a jerky shifting of the phase; 
the April 23, 1961 recording looks rather like a 
flight of five steps. It will be shown that each step 
is associated with a decrease in the altitude of a 
night reflection point on the ionosphere. Identical 
symbols were maintained for similar singularities 
on different steps, indexes referring to step numbers. 
These symbols are the same as those used previously 
in the recordings of the GBR signal and relate to 
the same singularities of the phase shift (sec. 7.2b). 

Reference is made here to figure 13, giving a 
schematic description of the model used. The 
ionospheric reflection points involved in these 
computations are to be found between the diagonal 



134 



lines in table 3 (sec. 6.2a). Times of sunrise and of 
the rising of the ionizing radiation as well as cor- 
responding zenithal distances of the sun were 
computed for each of the reflection points mentioned 
(table 13). The night reflection altitude had been 
selected at 84 km in section 6.2a. 

Table 14 is a systematic collection of the singu- 
larities noted in the recorded phase shifts of the 
NBA signal. They were assigned the letters H i} I u 
K iy L t , Oi, Mi (i= 1,2,3,4,5), i being the number of 
the step studied. 

The following considerations are based on table 14: 

(1) Contrary to the results obtained from the 
phase shift of the GBR signal, a first phase shift 
occurs at times H t prior to the shift at I t . Table 14 
shows that at the times H t the zenithal distances 
computed for the reflection points E\ vary between 
101 and 104° while the distances d from the reflection 
points E\ to the layer of ionizing radiation fluctuate 
between 60 and 150 km. On both recordings of 
September 8 and 12, I960, the zenithal distances 
calculated for the times B A (08.30 UT, 08.37 UT) 
had to be eliminated since the recordings are inter- 
rupted at 08.30 UT. The computation Avas made to 
show that these singularities belong to the third step 
(singularity 3 ) . 

(2) At times I t corresponding to the rising of the 
ionizing radiation at reflection point E on the re- 
cordings of the GBR signal (sec. 7.2b.), the zenithal 
distances related to point E\ vary between 99 and 
101°, the distances d to the layer formed by the 
ionizing radiations being between 6 and 93 km. 

(3) At the times t or M u the zenithal distances 
related to the corresponding day reflection points 



D 2 



T>\, D|, Dl, D? are near 90°. This result cor- 



roborates the findings on the recordings of the 
GBR signal. 

(4) The zenithal distances of the sun calculated 
at the intermediate times K± and L t are generally 
between 93 and 97°. 



(5) The energy A 2 of the NBA signal ((3) sec. 1) 
was shown below the April 23, 1961, recording as a 
function of time. For reasons of convenience, A 2 
was represented as the reciprocal function 1/A 2 due 
to the large variations of this function. This curve 
shows that: 

The decrease in the energy occurs only after 05 UT 
(time L{). For that time, the zenithal distance of 
the sun at the day reflection point D\ would be 
92°48', table 14. "The absorption curve of this 
initial step cannot be compared to the corresponding 
curve for Rugby, figure 5. 

Two relative minima are discernible at the times 
2 and M 2 . These refer to the second step. 

The third, fourth, and fifth minima occurring at 
the times 3 , 0* and is are associated to steps 
3, 4, 5. Zenithal distances computed with respect 
to the corresponding day reflection points fluctuate 
between 92 and 93°. 

The intermediate maxima in the energy occur 
generally between the times Hi and 7 Z ; these times 
precede the rising of the ionizing radiations of the 
sun at the night reflection point E[ only shortly. 

Recordings of September 7, 1960, Figure 6; September 11, 1960, Figure 
8; and April 22, 1961, Figure 10 at Sunset 

The shifting of the phase of the NBA signal is 
much smoother at sunset than at sunrise. The dif- 
ferent steps so clearly discernible at sunrise are 
smoothed out to such an extent that an analysis of 
the recordings becomes quite contingent. This is 
easily understood by referring to the phase recording 
of the GBR signal at sunset, figures 2 and 4. At- 
tention had already been drawn to the considerably 
longer time periods over which the phase shift took 
place at sunset as compared to sunrise. Applied to 
a multihop skywave, this means that the phase shift 
associated with tin 4 first reflection point at sunset is 
not terminated by the time the shifting associated 
with sunset at the second reflection point intervenes. 



Table 13 



Reflection 




Time of 


sunrise 




Time of incidence of the first ionizing sunray 


point or 
place 


8.9.60 
UT 


12.9.60 
UT 


23.4.61 
UT 


Z 


8.9.60 
UT 


12.9.60 
UT 


23.4.61 
UT 


Z 


Neuchatel 


5fc05"» 

4 57 


5 hHm 

5 03 


4 h 36"> 
4 24 


90°00' 
99 15 










E\ 


5'>07«i 


5 hl4m 


4 h 36 m 


97°33' 


D\ 


4 59 


5 04 


4 26 


98 27 


5 10 


5 16 


4 39 


96 33 


El 


6 28 


6 31 


5 59 


99 15 


6 38 


6 41 


6 10 


97 33 


D\ 


6 31 


6 36 


6 02 


98 27 


6 42 


6 47 


6 L3 


96 33 


E\ 


7 53 


7 56 


7 30 


99 15 


8 02 


8 05 


7 39 


97 33 


i>\ 


7 57 


8 00 


7 35 


98 27 


8 07 


8 10 


7 45 


96 L3 


E\ 


9 05 


9 07 


8 50 


99 15 


9 13 


9 15 


8. 58 


97 33 


D\ 


9 10 


9 12 


8 55 


98 27 


9 19 


9 21 


9 04 


96 33 


E\ 


10 08 


10 08 


9 59 


99 15 


10 15 


10 15 


10 06 


97. 33 


l\ 


10 13 


10 13 


10 05 


98 27 


10 21 


10 21 


10 13 


96 33 
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Table 14 





Sunrise 


Abscissa in 
the figures 


Reflection 
point 


8.9.60 


12.9.60 




23.4.61 




7, 9, and 10 


UT 


Z 


d 


UT 


Z 


d 


UT 


Z 


d 










km 






km 






km 


Hi 


E\ 


4^32"' 


103°13' 


115 


4 h 46m 


101°56' 


84 


4hl 4 m 


100°43' 


57 


h 


K\ 


4 58 


99 05 


25 


5 02 


99 25 


31 


4 29 


98 30 


15 


K x 




5 15 
5 29 


95 50 
93 12 




5 12 
5 39 


97 07 
92 36 


8 


4 46 

5 03 
5 26 


95 26 
92 48 
89 06 




Lx 






Mi 
























Ho 




5 57 

6 28 


104 31 
99 15 


150 

28 


6 00 
6 31 


104 35 

99 15 


152 

28 








h 


6 07 


97 57 


6 


K% 


D\ 
I)\ 
D\ 


6 44 
6 59 


96 11 
93 33 




6 38 


98 06 


22 


6 07 
6 19 
6 36 
6 55 


97 38 

95 40 
92 55 
89 33 


15 


L 2 






Oo 




7 13 

7 18 


91 54 
91 00 






M 2 
























m 


E] 


7 32 


103 18 


117 


7 29 


104 29 


149 


7 14 


102 09 


89 


h 


E\ 


7 47 


100 25 


50 


7 47 


101 01 


63 


7 33 


98 42 


18 


K 3 










8 08 


96 53 


4 








U 








7 52 

8 09 


95 11 
91 57 




o 3 


8 30 


91 58 




8 37 


91 24 














m 


E\ 


8 30 


106 50 


221 


8 37 


105 48 


188 


8 34 


102 36 


99 


h 


E\ 














8 43 


100 43 


57 


K* 


D\ 


Int 


err up ted rec 


3rd 


Int 


errupted record 


9 00 


97 24 


11 


o 4 


B\ 














9 21 


92 54 








Is 


E\ 














9 59 


99 15 


28 


K h 


^5 














10 12 


96 49 


4 


L- 
















10 29 

11 24 


92 51 
79 50 




Po 









An overriding of phase shifts occurs resulting in 
smoother curves for sunset than for sunrise. 

Nevertheless, a trial was attempted. Although 
results obtained may not be very well established, 
there are indications that the schematic model de- 
picted in figure 14 may be fairly representative of 
reality. The reflection points considered are to be 
found between the diagonal lines of table 6. 

Table 15 contains the times of rising of the ionizing 
radiations of the sun at each of the ionospheric re- 
flection points mentioned (see fig. 14 and table 6) 
together with sunrise times at these points and re- 
lated zenithal distances. The altitudes of the night 
reflection points had been selected at 84 km (sec. 
6.2a). 

The times designated by the symbols A% A iy B u 
Ct, F t , F% and G* (table 16) were selected whenever 
phase shift singularities were more or less apparent. 
Not all were chosen due to the fairly even variation 
of the recorded curves. 



Examination of table 16 and the recordings, figures 
6, 8, and 10 leads to the following considerations: 

(1) Shifting of the phase of the NBA signal usually 
begins at times A t (i= 1,2,4) with zenithal distances 
of the sun with respect to the reflection points Df 
near 90°. 

(2) On the April 22, 1961 recording, figure 10, the 
initial phase shift occurs earlier at time A*, when the 
zenithal distance of the sun is only 85,5°. 

(3) The singularities found at times Bi and B 2 are 
peculiar to the recording of September 11, 1960, 
figure 8. Relevant zenithal distances are near 94°. 

(4) At the times d (i= 1,2,3,5) zenithal distances 
fluctuate between 97 and 100°, distances d from the 
night reflection point to the layer formed by the 
ionizing radiations between and 29 km. It can 
therefore be assumed that the times d are approxi- 
mately the setting times of the ionizing radiation at 
the night reflection points E\< respectively El, E$, 
and m. 
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Table 15 



Reflection 


'lime of sunset 


Time of incidence of the last ionizing sunray 


point or 

place 


7.9.60 

UT 


11.9.60 
UT 


22.4.61 
UT 


Z 


7.9.60 
UT 


Ll.9.60 

UT 


22.4.61 
UT 


Z 


Neuchatel 


17^56™ 
19 34 


17147m 
19 26 


18h25"> 

20 06 


90°00' 

98 27 










D\ 


19'>22m 


19hl4m 


19^53™ 


96°33' 


E\ 


19 42 


19 33 


20 14 


99 15 


19 32 


19 23 


20 03 


97 33 


D\ 


21 02 


20 54 


21 30 


98 27 


20 51 


20 43 


21 19 


96 33 


E\ 


21 07 


21 00 


21 36 


99 15 


20 58 


20 50 


21 25 


97 33 


D\ 


22 07 


22 01 


22 28 


98 27 


21 58 


21 51 


22 18 


96 33 


n 


22 11 


22 04 


22 33 


99 15 


22 02 


21 56 


22 24 


97 33 


D\ 


22 56 


22 52 


23 11 


98 27 


22 48 


22 43 


23 02 


96 :y,\ 


E* 


23 00 


22 55 


23 14 


99 15 


22 52 


22 47 


23 06 


\)7 :ya 


D\ 


23 37 


23 33 


23 45 


98 27 


23 29 


23 25 


23 37 


96 33 


F-l 


23 39 


23 35 


23 48 


99 15 


23 32 


23 28 


23 40 


97 33 



Table 16 



Sunset 



Abscissa in 


Reflection 
point 


7.9.60 


11.9.60 


22.4.61 


the figures 
6, 8, and 10 


UT 


7. 


(1 


UT 


Z 


(I 


UT 


Z 


d 


Ax* 


E\ 

E\ 

K 

E\ 
El 






km 






km 


18h38m 


85°36' 


km 


Ax 


l 8 h4 7 m 


90°39' 












B x 




19M)0 m 
19 34 


94°08' 
99 25 










Ci 


19 27 

20 14 
20 14 


96 48 
104 20 
89 48 




31 


20 18 


99 51 


39 


Fx 


145 




A<> 














B 2 




20 32 
20 56 


94 31 
96 49 










C 2 


21 02 

22 02 


98 22 
97 33 


13 











C-i 




22 31 


98 52 


21 


Fi 


22 24 
22 24 


103 14 
92 11 


115 




A\ 








22 31 

23 49 


89 46 
99 29 




C 5 










32 


Fh 


23 45 

24 06 


100 41 
105 41 


56 
184 


23 46 


101 53 


82 




F$* 








G 5 * 








25 21 
25 40 


120 23 
124 31 


830 


Go 














1080 



















(5) At the times Fi, F 3 , and F 5 but only on the Sep- 
tember 7, figure 6 and September 11, I960, figure 8 
recordings zenithal distances with respect to the 
night reflection points El, respectively E%, E%, vary 
between 101 and 104°; the distances d would then 
be of the order of 100 km. 

(6) The times Ft figure 6, G\ and G h figure 10 
were indicated and relevant zenithal distances and 
distances d were computed, however, without 
permitting any significant conclusion. 

(7) The energy of the NBA signal ((3) sec. 1), 



was plotted on the April 22, 1961 recording as a 
function of time (see sec. 7.3a, item 5). The 
energy decreases after about 17 UT when the 
zenithal distance of the sun with respect to the day 
reflection point D\ becomes smaller than 85°. At 
time At (Z=85°36 / ), the function 1/Af diminishes 
rapidly. A weak maximum is discernible at time 
C\ as well as a bend at time Aj = 3 . This absorption 
curve shows a diminished energy during the time 
period between sunset at reflection point D\ and 
sunset at the last reflection point D\. 
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8. Conclusions 

Computation of angles of incidence <j> on the ground 
as well as the deduction of the altitude of the night 
reflection points were made possible by determination 
of the propagation modes of the GBR signal (Rugby) 
and of the NBA signal (Balboa). The altitude of 
the day reflection points was taken from literature. 
Angles of incidence on the ground for the NBA 
signal vary between 0°27' and 1°14'. These angles 
are much smaller than those derived from actually 
accepted mathematical theories [Wait 1961, Wait 
and Conda 1961]. 

The detailed study of these propagation modes 
was made possible by introducing the zenithal 
distances of the sun. A good accord was shown to 
exist between the zenithal distances computed for 
times selected on recordings with respect to the 
reflection points determined by the model, figures 
13 and 14. The validity of the model is corrobo- 
rated by these relations. Zenithal distances of the 
sun show a recurrence of these values on various 
recordings: 

(1) Z^90° is the zenithal distance at which the 
altitude of the reflection point begins to increase at 
sunset. The altitude of the reflection point reaches 
the day value at sunrise. 

(2) Z^97° is the zenithal distance at which the 
onset of the ionizing radiation occurs at the night 
reflection points. An increased phase shift is 
frequently observed at this time. 

(3) Z^103° (for sunrise on the Balboa recordings 
only) is the zenithal distance at which the layer 
formed by the ionizing radiations is still about 100 
km distant from the night reflection point. This 
time marks the beginning of the phase shift of the 
NBA signal and is identical with the beginning 
decrease of the night reflection point altitude. The 
decrease of the altitude of the night reflection pre- 
cedes the onset of the ionizing radiation of the sun 
and substantiates previous statements [Rieker 1960, 
sec. 3.5.2.3.]. 

The fact that a similar phenomenon does not occur 
on the GBR signal transmission is probably due to 
to the higher altitude of the night reflection point 
(see sec. 4.4). 

At sunrise, the behavior of the energy generally 
reveals a pronounced absorption at the time at 
which the phase shift, resulting from the decrease 
of the altitude of one ionosphere reflection point is 
about half completed. This absorption maximum 
occurs for zenithal distances of the sun varying 
between 92 and 95°, figures 5 and 10. 

At sunset, the energy reveals a minor maximum 
around 17 UT, figure 4, for a zenithal distance near 
94°. The NBA signal recorded in figure 10 merely 
shows a weakening of the energy during sunset on 
the path Neuchatel-Balboa. 



The absorption curves plotted under the record- 
ings of figures 4,5,10, and the zenithal distances 
included in the tables are only of a descriptive nature. 
The study of their behavior has facilitated the under- i 
standing of the relation existing between sunrise 
and sunset in the ionosphere and the propagation of 
longwaves. Now that the general behavior of 
this relation appears to be known qualitatively and 
that the close relationship with the energy and the I 
phase shift has become obvious, a study of the de- 
tailed quantitative behavior of the energy would be 
very interesting. This study would require, so 
it seems, a more direct method of absorption 
measurement. 



We record our gratitude to Dr. J. Bonanomi, Direc- 
tor of the Neuchatel Cantonal Observatory, who 
made available the recordings as well as relevant 
technical data. 
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